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Abstract. We study for what specific values of the theoretical parameters the axion can form 
the totality of cold dark matter. We examine the allowed axion parameter region in the light of 
recent data collected by the WMAP5 mission plus baryon acoustic oscillations and supernovae 
[1] , and assume an inflationary scenario and standard cosmology. We also upgrade the treatment 
of anharmonicities in the axion potential, which we find important in certain cases. If the Peccei- 
Quinn symmetry is restored after inflation, we recover the usual relation between axion mass 
and density, so that an axion mass m a = (85 ± 3) /xeV makes the axion 100% of the cold dark 
matter. If the Peccei-Quinn symmetry is broken during inflation, the axion can instead be 100% 
of the cold dark matter for m a < 15 meV provided a specific value of the initial misalignment 
angle 8i is chosen in correspondence to a given value of its mass m a ■ Large values of the Peccei- 
Quinn symmetry breaking scale correspond to small, perhaps uncomfortably small, values of 
the initial misalignment angle 8i. 



1. Introduction 

About 84% of the non-relativistic matter in the Universe is in the form of cold dark matter 
(CDM)[1], whose nature is still to be discovered. One of the most promising hypothetical 
candidates that could account for the CDM observed is the axion [2, 3]. 

The properties of the axion as the CDM particle have been studied in various papers [4, 5]. 
We examine the possibility that the axion accounts for 100% CDM in the light of the WMAP5 
mission, baryon acoustic oscillations (BAO) and supernovae (SN) data. We also upgrade the 
treatment of anharmonicities in the axion potential, which we find important in a specific 
region of the parameter space. The axion parameter space is described by three quantities, 
the PQ energy scale /„, the Hubble parameter at the end of inflation Hj and the axion initial 
misalignment angle Q{. 



2. Properties of the axion 

The solution to the strong-CP problem by Peccei and Quinn [3] introduces a new U(1)pq 
symmetry in the theory of strong interactions that spontaneously breaks at a temperature of 
order the PQ energy scale, T ~ f a . The axion is the pseudo Nambu-Goldstone boson originating 
from the breaking of U(1)pq. We denote by a(x) the axion field and 6{x) = a(x)/f a the 
misalignment field. At a temperature of order the QCD transition T ~ h-QCD an effective 
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potential V{9) arises through instanton effects [6] given by 



V(9)=ml{T)fl(l-co S 6). (1) 
The temperature-dependent axion mass m a (T) is given by 

ma{T)= U a 6(A) 4 T > Aqcd, (2) 
\m a T < Aqcd- 

We set Aqcd = 200MeV, b = 0.018 [5]. Setting f a>l2 = / a /10 12 GeV, the axion mass at zero 
temperature is 

m a = 6.2^eVf-l 2 /N. (3) 
The integer TV" represents the U(1)pq color anomaly index; here N = 1. 

3. Axion and cosmology 

Setting Tqh = Hi/2tt, axion cosmology is different if the PQ symmetry breaks after inflation, 
f a < Tch (Scenario I), or during inflation, f a > Tch (Scenario II). Over a Hubble volume, 9i 
takes different values in Scenario I, while in Scenario II 9i takes only one random value. 

The subsequent evolution of the zero mode of the dynamical field 9{x) in a flat Friedmann- 
Robertson- Walker metric is given by 

+ 3H{T) + l_^l = . (4) 

Here a dot indicates a derivative with respect to time and H(T) is the Hubble expansion rate. 
For T 3> Aqcd the misalignment field is frozen at a constant value, 9 = 9i = const, where 6i 
is the initial misalignment angle. For m a (T) ~ 3H(T) the axion field begins to oscillate in the 
V{9) potential. The axion energy density-to-critical density is computed as 

n mis h 2 = f 0.236<^m)> f a % fa < fa, (g) 

(0.0051(^/(0.))/^, /a >/ 0j 

with f a = 9.91 x 10 16 GeV. The function f(9i) accounts for anharmonicities in the axion potential, 
i.e., deviations from a pure harmonic potential when 9\ is no longer smaller than one. We look 
for an analytic function f(0i) that interpolates between previous numerical results [7]. We take 
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(6) 



The function f(9i) in eq.(6) shows the correct asymptotic behaviors for 9i — > ir and 9i — > 0. 



4. Results 

Figure 1 shows the region of the axion parameter space in which the axion is 100% CDM, or 
ft a = ^cdm, where VLcdm is the CDM density in units of the critical density. The PQ scale f a 
is bounded from below by measurements on white dwarfs cooling time, /„>4x 10 8 GeV, while 
the Hubble rate at the end of inflation Hi is bounded from above by non-detection of tensor 
modes, tfj < 6.3 x 10 14 GeV. For f a < T GH (Scenario I), the axion can be 100% CDM only if 
f a = (7.27 ± 0.25) x 10 10 GeV. The function jf(0j) in eq.(6) affects the results by increasing the 
value of (9ff(9i)). In Scenario II, f a > Tqh, the parameter space on the top left is bounded by 



non-detection on axion isocurvature fluctuations that yield f a < 4.17 x 10~ 5 Hj Q{. The upper 
kink in the leftmost boundary is due to the change in the dependence of the axion mass with 
temperature, eq.(2), while the lower kink is due to the presence of anharmonicities, eq. (6). In 
the allowed region of Scenario II, for a given f a the axion is 100% CDM only if a correspondent 
value 6{ = 0i(f a ) is selected. This value of 6i does not depend on Hi. We show the function 
0i = Oi(fa) in Figure 2, where we see that a value 9i < 6i(f a ) for a given f a makes the axion 
population sub-dominant compared to the CDM observed, while a value 9i > 0i(f a ) is over 
abundant. We included the region excluded by ADMX1 axion direct search and future reaches 
by the Planck satellite and future cavity searches. 
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Figure 1. Region of axion parameter space where 
the axion is 100% of the cold dark matter. The axion 
mass scale on the right corresponds to eq. (3) with 
N = 1. When the PQ symmetry breaks after inflation 
{fa < Hi/2tt), the axion is 100%CDM if f a = (7.27 ± 
0.25) x 10 10 GeV, or m a = (85 ± 3) ^eV, which is 
the narrow horizontal window shown on the right. If 
the axion is present during inflation (f a > Hi/2rr), 
axion isocurvature perturbations constrain the parameter 
space to the region on the top left, which is marked 
by the values of Qi necessary to obtain 100% of the 
CDM density. Other bounds indicated in the figure come 
from observations of white dwarfs cooling times and the 
non- observation of tensor modes in the CMB spectrum. 
Dashed lines and arrows indicate the future reach of 
the PLANCK satellite and the ADMX and CARRACK 
cavity searches. 
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Figure 2. The misalignment angle 8i necessary 
for the axion to be 100% of the cold dark matter 
in Scenario II (f a > Hi [lit), as a function of f a . 
Above the curve, fi a > Qcdm- For f a > 10 17 GeV, 
ft w 0.001(/ a /10 17 GeV) _3/4 . 



